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Abstract 

We look at geometrical and arithmetical patterns created from a finite 
subset of Z n by diffracting waves and bipartite graphs. We hope that this 
can make a link between Motives and the Melting Crystals/Dimer models 
in String Theory. 

Introduction. 

Why is it that, occasionally, mathematicians studying Motives and physicists 
searching for a Theory of Everything seem to be looking at the same examples, 
just from different angles ? Should the Theory of Everything include properties 
of Numbers? Does Physics yield realizations of Motives which have not been 
considered before in the cohomological set-up of motivic theory? 

Calabi-Yau varieties of dimensions 1 and 2, being elliptic curves and Re- 
surfaces, have a long and rich history in number theory and geometry. Calabi- 
Yau varieties of dimension 3 have played an important role in many develop- 
ments in String Theory. The discovery of Mirror Symmetry attracted the at- 
tention of physicists and mathematicians to Calabi-Yau's near the large complex 
structure limit [131 12U| . Some analogies between String Theory and Arithmetic 
Algebraic Geometry near this limit were discussed in |16l 1171 IT%| . Recently new 
models appeared, called Melting Crystals and Dimers |14l II 1| . which led to in- 
teresting new insights in String Theory, without going near the large complex 
structure limit. The present paper is an attempt to find motivic aspects of these 
new models. We look at geometrical and enumerative patterns associated with 
a finite subset 21 of Z n . The geometry comes from waves diffracting on 21 and 
from a periodic weighted bipartite graph generated by 21. The latter is related 
to the dimers (although here we can not say more about this relation). Since 
the tori involved in these models are naturally dual to each other there seems 
to be some sort of mirror symmetry between the diffraction and the graph pic- 
tures. The enumerative patterns count lattice points on the diffraction pattern, 
points on varieties over finite fields and paths on the graph. They are expressed 
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Figure 1: Diffraction pattern for 21 = {(1,0), (0,1), (-1,-1)} c Z 2 , allc 3 = l. 

through a sequence of polynomials Bat(z) with coefficients in Z and via a limit 
for zeC: 

Q(z)= lhn |B„(z)|- w -\ (1) 

Limit formulas like Q appear frequently and in very diverse contexts in the 
literature, e.g. for entropy in algebraic dynamical systems in [7] Theorem 4.9, 
for partition function per fundamental domain in dimer models in |11| Theorem 
3.5, for integrated density of states in [2] p. 206. Moreover, Q(z) appears as 
Mahler measure in , as the exponential of a period in Deligne cohomology in 
|6l 115) . and in instanton counts in |18| : see the remark at the end of Section 

With some additional restrictions 21 provides the toric data for a family of 
Calabi-Yau varieties and various well-known results about Calabi-Yau varieties 
near the large complex structure limit can be derived from the Taylor series 
expansion of logQ(z) near z = oo; see the Remark at the end of Section [S] In 
the present paper we are not so much interested in the large complex structure 
limit. Instead we focus on the polynomials Bn(z) and the limit formula Q). 
This does not require conditions of 'Calabi-Yau type'. 

When waves are diffracted at some finite set 21 of points in a plane, the 
diffraction pattern observed in a plane at large distance is, according to the 
Frauenhofer model, the absolute value squared of the Fourier transform of 21. 
There is no mathematical reason to restrict this model to dimension 2. Also 
the points may have weights > 1. So, we take a finite subset 21 of Z" and 
positive integers c a (a G 21). These data can be summarized as a distribution 
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V = J^aea c a^a ) wnere S a denotes the Dirac delta distribution, evaluating test 
functions at the point a. The Fourier transform of T> is the function T>(\) = 
Saesi c a e~ 27 ™^ a ^ on 1"; here (,) is the standard inner product on 1". The 
diffraction pattern consists of the level sets of the function 

\V(t)\ 2 = Cac b e 27 " (t < a - b> = J2 c a c b cos(2^(t,a-b)). 
a,bea a,bea 

This function is periodic with period lattice A v dual to the lattice A spanned 
over Z by the differences a — b with a, b € 21. Throughout this note we assume 
that A has rank n. Looking at the intersections of the diffraction pattern with 
the lattices ^A v we introduce the enumerative data 

multjvM := tt{t G ^A v /A v | \V(t)\ 2 = r } for N G N, r G K , (2) 
and use these to define polynomials B^(z) as follows: 
Definition 1 

Bjv(2):=n(*- r ) mUlMr) - ( 3 ) 

rGR 

One could also introduce the generating function F(z,T) := X^jvgn &n(z)T n , 
but except for the classical number theory of the case n — 1 (see Section , 
and the observation that Formula gives Q(z) as the radius of convergence 
of F(z,T) as a complex power series in T, we do not yet have appealing results 
about F(z,T). 

For the graph model we start from the same data: the finite set 21 C Z", 
the weights c a (a G 21) and the lattice A spanned by the differences a — b with 
a, b G 21. We must now assume that 

21 n A = 0. 

One can then construct a weighted bipartite graph T as follows. Bipartite 
graphs have two kinds of vertices, often called black and white. The set of black 
vertices of T is A. The set of white vertices of T is 21 + A. Note that 21 + A is just 
one single coset of A in Z™ . In F there is an (oriented) edge from vertex Vi to 
vertex V2 if and only if Vi is black, V2 is white and V2 — Vi G 21. If V2 — vj = a G 21 
the edge is said to be of type a and gets weight c 3 . 

The graph T is A-periodic and for every N G N one has the finite graph 
Tjv := r/ 7VA, which is naturally embedded in the torus R n / NA- By a closed 
path of length 2k on T or we mean a sequence of edges (ei, e%, . . . , e2fc_i, e2fc) 
such that for i = 1, . . . , k the intersection e2i_i He2i contains a white vertex and 
e2i H e2i+i contains a black vertex; here e2fc+i = ei. By the weight of such a 
path we mean the product of the weights of the edges ei, ea, • . . , e2fc-i, £2k- We 
denote the set of closed paths of length 2k on Tn by rjy(2fc). Enumerating the 
closed paths on Tat according to length and weight we prove in Section [3] that 
this leads to a new interpretation of the polynomials Bn(z): 
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Figure 2: A fundamental parallelogram of the lattice 3 A and a piece of the 
bipartite graph T for 21 = {(1,0), (0,1), (-1,-1)} C I?, and all c 3 = 1. By 
identifying opposite sides of the parallelogram one obtains the graph IV 



Theorem 1 

B N (z) = z Nn exp 1-^2 E weight( 7 )z- fc ) . (4) 

\ fc>l76r„(2fc) / 



Formulas © and transfer the enumerative data between the two models. 
We pass to algebraic geometry with the Laurent polynomial 

W( Xl ,...,x n ) = Y, c a c b x a - b &1[xf x ,...,xt l ], (5) 

a,b£2l 

which satisfies |2?(t)| 2 = W{e 2mt \ . . . , e 27 ™*"); here x A := if A = 

(Ai, . . . , A„) S A. For A" S N let /^at denote the group of AT-th roots of unity 
and let /x^ := Hom(A,/i/v) be the group of homomorphisms from the lattice A 
to /ipf. Thus the defining formula can be rewritten as: 

Btt(z)= LI s - W(x)) . (6) 

Written in the form ijBJ the polynomials Bn(z) appear as direct generalizations 
of quantities introduced by Lehmer |12j for a 1-variable (i.e. n = 1) polynomial 
W(x). Using JJjJ one can easily show (Proposition^ that the polynomials Bn(z) 
have integer coefficients and that Bn'(z) divides Bn(z) in Z[z] if N' divides A~ in 
Z. Thus for h € Z also Bjv(^) is an integer. Lehmer was particularly interested 
in the prime factorization of these integers in case n = 1 |12l El- Also 
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for general n > 1 these prime factorizations must be interesting, for instance 
because they relate to counting points on varieties over finite fields; see Section 
0| for details. Thus prime factorization gives a third occurrence of Bn(z) in 
enumerative problems, related to counting points on varieties over finite fields. 
Q{z) appears in ^0^1 as Mahler measure with ties to special values of L- 
functions. It would be nice if the limit formula (|TJ together with the prime 
factorization of the numbers Bjv(z) (with z £ Z) could shed new light on these 
very intriguing ties. 

In Section Q] we study the density distribution of the level sets in the diffrac- 
tion pattern. Passing from measures to complex functions with the Hilbert 
transform we find one interpretation of Q(z), Bjv(z) and QJ. In Section\^we 
briefly discuss another interpretation in connection with the spectrum of a dis- 
cretized Laplace operator. In Section\^we prove Theorem^ In Section ^] we 
pass to toric geometry, where the diffraction pattern reappears as the intersec- 
tion of a real torus with a family of hypersurfaces in a complex torus and where 
logQ(z) becomes a period integral, while the prime factorization of Bn(z) for 
z G Z somehow relates to counting points on those hypersurfaces over finite 
fields. In Section [5| we discuss sequences of integers which appear as moments 
of measures, path counts on graphs and coefficients in Taylor expansions of so- 
lutions of Picard-Fuchs differential equations. Finally, in Section we present 
some concrete examples. 

1 The diffraction pattern. 

The function |£>(t)| 2 is periodic with period lattice A v dual to the lattice A: 

A v := {t e R" | (t, a - b) e Z , Va, b e 21} , 
A := Z-Span{a - b | a, b G 21}. 

Throughout this note we assume that the lattices A and A v have rank n. 

Because of this periodicity |£>(t)| 2 descends to a function on R" / A v . Defin- 
ing for t S R" the function e t : R" — > C by e t (v) = e 2,ri < t ' v ) we obtain an 
isomorphism of real tori 

R«/ A v ~ U A , t^e t , (7) 

where U A := Hom(A,U) is the torus of group homomorphisms from the lattice 
A to the unit circle U:={a;eC||x| = l}. Recall that a group homomorphism 
ip : A — > U induces an algebra homomorphism -0* from the group algebra C[A] 
to C. Thus C[A] is the natural algebra of functions on U A and -0* evaluates 
functions at the point ip of U A . The inclusion A C Z™ identifies C[A] with 
the subalgebra of the algebra of Laurent polynomials Cfof 1 , . . . , x^ 1 ], which 
consists of C-linear combinations of the monomials x A := Oi=i w ith ^ = 
(Ai, . . . , A„) £ A. Thus, via J7J), the function |2?(t)| 2 coincides with the Laurent 
polynomial W(xi, . . . , x n ) defined in JSJ). 
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Positivity of the coefficients c a implies that the function |D(t)| 2 attains its 
maximum exactly at the points t G A v . In terms of the torus U A and the 
function W : U A — > R this means that W attains its maximum exactly at the 
origin 1 of the torus group U A : 

Vx 6 U A \ {1} : W(x) < W(l) = C 2 with C := ^ c a . 

aea 

Some important aspects of the density distribution in the diffraction pattern 
are captured by the function 

We view the derivative dV(r) of V as a measure on R. In our analysis it will 
be important that the measure dV(r) is also the push forward of the standard 
measure dt\ dt^ ■ ■ ■ dt n on R™ by the function |2?(t)| 2 . 

Another insight into the diffraction pattern comes from its intersection with 
the torsion subgroup of U A . For N G N let /ijv C U denote the group of iV-th 
roots of unity. Then the group of TV-torsion points in U A is Hjj := Hom(A, /ijv) 
and (J2J) can be rewritten as 

multjv(r) = |j(/4 n W- 1 ^)) for N G N, r G R . 

Moreover we set, in analogy with JSJ, 

Mr) := T^tt{x G a4 I W(x) <r} for N G N, r G R. 

iv " 

The derivative of the step function V/v : K — » R is the distribution 

dT/jv (r) = iV~ n ^ mult w (r) 5 r , (9) 

r 

which assigns to a continuous function / on 1R the value 

[ f(r)dV N (r) := N~ n ^ multjv(r) f(r) = N~ n £ /(W(x)) . (10) 

One thus finds a limit of distributions 

lim dV N (r) = dV(r); (11) 

TV— »oo 

by definition, this means that for every continuous function / on R 

lim / f(r)dV N (r) = [ f(r)dV(r). (12) 
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Measure theory is connected with complex function theory by the Hilbert 
transform. The Hilbert transform of the measure dV(r) is the function —-H(z) 
defined by 

H(z) := / — !— dV{r) for zeC\I; (13) 



here J := {r e M | < V(r) < 1} C [0, C 2 ] is the support of the measure dV(r). 
The measure can be recovered from its Hilbert transform because for tq S K 

^ (ro) = i £ eteo (H(ro + Je) - H(r °- ?e)) - 
Another way of writing the connection between dV(r) and H(z) is 

^ ^ /(z)H(*)<fe = j^(±-J |^<fe) dV(r) = jf /(r)dF(r) (14) 

for holomorphic functions / defined on some open neighborhood U of the interval 
I in C and closed paths 7 in f \ I encircling I once counter clockwise. 
Next we consider the function 

Q : C \ 1 — ► M >0 , Q(z) := exp f- jf log |z - r|dF(r)^ . (15) 

This function satisfies 4^ logQ(z) = — H(z) and thus ifTH can be rewritten as 

^-.j f{z)d\ogQ{z) = Jj(r)dV(r). 

This means that, at least intuitively, the functions Q(z) and e^ v ^ correspond 
to each other via some kind of comparison isomorphism. 

In order to find the analogue of 1)11(1 in terms of functions on C \ X we apply 
(fTUf) to the function f(r) = log \z - r\ on K with fixed z e C \ 1: 

[ log\z~r\dV N (r) = N- n \og\B N (z)\ (16) 

where B N {z) = l\ r€l (z - r ) mult «M = U^( z ~ W ( x )) as in © and El- 
Combining i|12[) . (|15|l and (|16|) we find the limit announced in QJ: 

Proposition 1 Q(z) = lim |B7v(;z)| _Ar for every z G C\X . ■ 

AT— *oc 



7 



2 The Laplacian perspective. 

Convolution with the distribution T> gives the operator D/(v) := E] c a /(v — a) 

aiE2l 

on the space of C- valued functions on W 1 . Let D/(v) := E] c a /(v + a) and 
A := DD , A/(v) = E c a c b /(v + a - b) . 

a,b<E2l 

For a sufficiently differentiable function / on R™ the Taylor expansion 
A/(v) = C 2 /(v) + | E E c ^( a * ' - b ^d^ {y) + ' ' ' 

shows that the difference operator A — C 2 is a discrete approximation of the 
Laplace operator corresponding to the Hessian of the function |2?(t)| 2 at its 
maximum. 

Remark. In |5] Gieseker, Knorrer and Trubowitz investigate Schrodinger equa- 
tions in solid state physics via a discrete approximation of the Laplacian. In 
their situation the Schrodinger operator is the discretized Laplacian plus a pe- 
riodic potential function. So from the perspective of [H] the present note deals 
only with the (simple) case of zero potential. On the other hand we consider 
more general discretization schemes and possibly higher dimensions. 

We now turn to the spectrum of A. For t £ K™ the function e t : M™ — * C 
given by e t (v) = e 27 ™^'^ is an eigenfunction for A with eigenvalue |2?(t)| 2 : 

Ae t (v) = E c aCbe 2 "^ v + a - b > = p(t)| 2 e t (v). 
a.bea 

Take a positive integer N. The space of C-valued C°°-functions on M. n which 
are periodic for the sublattice NA of A is spanned by the functions e t with t in 
the dual lattice ^A v . The characteristic polynomial of the restriction of A to 
this space is therefore (see © and ©) 

n (z-\v(t)\ 2 ) = n^-^w) = n(*- r ) muit " (r) = b ^)- 

te^A v /A v xe "« rm 

With 0) and Ml)) the measure dV(r) can now be interpreted as the density 
of the eigenvalues of A on the space of C-valued C°° -functions on I" which are 
periodic for some sublattice NA of A. 
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3 Enumeration of paths on a periodic weighted 
bipartite graph; proof of Theorem [U. 

In this section we prove Theorem ^ Recall from the Introduction just before 
Theorem n the various ingredients: the finite set 21 C the weights c a , the 
lattice A and the graphs T and Tjv- Recall also the closed paths on Tjv, their 
lengths and weights, and the set rjv(2/c) of closed paths of length 2k on T^. 

Consider a path (ei, e2, . . . , e2k-i, &ik) on T with edge going from black to 
white if i is odd, respectively from white to black if i is even. Let s denote the 
starting point of the path (i.e. the black vertex of edge ei). Let for j = 1, . . . , k 
edge e2j-i be of type 3j and edge e-ij of type bj. Then the end point of the 
path (i.e. the black vertex of e^k) is s + 2j=i ( a j ~ bj). The weight of the path 
is rX/_i c 3j c b r The path closes on Tn if and only if X^=i( a j — bj) 

Next recall from © that IF(x) = J2 a bea c a CbX a ~ b and set 



m W : = N- n W ( K ) k ■ 



xG/V^, 

So ml is the sum of the coefficients of those monomials in W(x) k with exponent 
in 7VA. In view of the above considerations mi is therefore equal to the sum 
of the weights of the paths on T which start at s, have length 2k and close in 
Tjv. Since on Tjv there are N n black vertices and on a path of length 2k there 
are k black vertices we conclude 



= weight( 7 ). (18) 

7Grjv(2fe) 



From © one sees for \z\ > C 
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m (iV) 



N- n \ogB N (z) = logz + N~ n J2 log(l- ^(x)z- 1 ) = logz-J2 , 

(19) 

Combining (jl%|) and (fTTijl we find 



B N {z) = z N " exp [ weight^) 

fe>i7erjv(2fe) 



z~ k 



This finishes the proof of Theorem ^ B 

Remark. In the Laplacian perspective A/^m^ is the trace of the operator 
A k on the space of C-valued C°°-functions on M. n which are periodic for the 
sublattice ATA of A. The polynomial B^(z) is the characteristic polynomial 
of A on this space. Formula l|19fl gives the well-known relation between the 
characteristic polynomial of an operator and the traces of its powers. 
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Remark. One may refine the above enumerations by keeping track of the 
homology class to which the closed path belongs. That means that instead 
of (|17|l one extracts from the polynomial W(x) the subpolynomial consisting 
of terms with exponent in NA. Such a refinement of the enumerations with 
homology data appears also in the theory of dimer models (cf. H2]), but its 
meaning for the diffraction pattern is not clear. 

4 Algebraic geometry. 

The polynomial Bjy(z) = Ilxe^ A ( z ~ ^( x )) nas coefficients in the ring of integers 
of the cyclotomic field Q(/ijv) and is clearly invariant under the Galois group 
of Q(/ijv) over Q. Consequently, the coefficients of Bn(z) lie in Z. The same 
argument applies to the polynomial Bn(z)Bn'(z)^ 1 = Ylx&n A \n A ( z ~ ^( x )) ^ 
N' divides N. Thus we have proved 

Proposition 2 For every N € N the coefficients of B^(z) lie in Z. If N' 
divides N in Z, then Bn'(z) divides Bjv(z) in 1\z\. M 

Fix a prime number p and a positive integer v G Z>o- Let W(F p ^) denote 
the ring of Witt vectors of the finite field W p * (see e.g. [3]). So, W(F p v) is a 
complete discrete valuation ring with maximal ideal pW(F p v) and residue field 
F p ^. The Teichmiiller lifting is a map r : F p « — > W(F p ^) such that 

x = t(x) mod p , i~(xy) — t(x)t(ii) \/x, y € F p k . 

Every non-zero x € F p >^ satisfies 

x^- 1 = 1. 

Thus there is an isomorphism /i^-i ~ F*„. Such an isomorphism composed 
with the Teichmiiller lifting gives an embedding j : fi p »-i W(F P *). Thus for 
x e we get 

W(j(x))eW(W p »). 
Recall the p-adic valuation on Z: for k E Z, k ^ 0: 

v p (k) := max{w £ Z | p v divides k} . 

Proposition 3 For p, v as above and for z € Z f/ie p-adic valuation of the 
integer Bp^_i(z) satisfies 

v p {B r ^(z)) > m G W I W(0 = 2 *nF p „ } . (20) 
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Proof: From © we obtain the product decomposition, with factors in W(F p ^ ) , 
The result of the proposition now follows because 



Remark. In l|3UI) we give an example showing that in i|2U[l we may have a strict 
inequality. 

Remark about the relation with Mahler measure and L-functions. 

The logarithmic Mahler measure m(i 7 ') and the Mahler measure M(F) of a Lau- 
rent polynomial F(x±, . . . , x n ) with complex coefficients are: 

l y y dx\ dx 
m i F ) : = ln - Vt ff log|F(n,...,af 7l )| 

M(F) := exp(m( J F)) . 

Boyd U gives a survey of many (two-variable) Laurent polynomials for which 
m(i 7 ') equals (numerically to many decimal places) a 'simple' non-zero rational 
number times the derivative at of the L-function of the projective plane curve 
Zp defined by the vanishing of F: 

m(F)-Q* = L'(Z F ,0)-Q*. (21) 

Deninger |Hj and Rodriguez Villegas JH] showed that the experimentally ob- 
served relations l|21JI agree with predictions from the Bloch-Beilinson conjec- 
tures. Rodriguez Villegas [R^ provided actual proofs for a few special examples. 

Since the measure dV(r) is the push forward of the measure dt\ dt 2 . . . dt n 
on E" by the function |2?(t)| 2 , one can rewrite Formula (|15fl as: 

-logQ(z) = -K- I log|*-inzi,---,Zn)| — (22) 
(2m) n J Vn xi x 2 x n 

On the right hand side of H22J) we now recognize the logarithmic Mahler measure 
of the Laurent polynomial z — W(x%, . . . , x n ) € C[xf 1 , . . . , a;^ 1 ]. 

For fixed z G Z Formulas and l|2(J|) provide a link between Q(z) and 
counting points over finite fields on the variety with equation W{x\, . . . , x n ) = z. 
It may be an interesting challenge to further extend these ideas to a proof of a 
result like (|2*TJ) . 



11 



5 Moments. 

Important invariants of the measure dV(r) are its moments (k £ Z>o): 

m fe := J R r k dV(r) = /* . . . £ \V(t x , . . .., t n )\ 2k dh . . . dt n 

= constant term of Fourier series \V(t\, . . . ,t n )\ 2k (23) 
— constant term of Laurent polynomial W(x%, . . . , x n ) k . 

The relation between the moments and the functions H(z), Q(z) defined in l|13[l 
and (0 is: for z £ R, z > C 2 , 



H(z) = Y, m kz- k ~\ Q(z) = z^exp K]^ z - fc | . (24) 

k>0 \k>l 



It is clear that the moments of dV(r) are non-negative integers. They 
satisfy all kinds of arithmetical relations. There are, for instance, recurrences 
like 1)28(1 and congruences like the following 

Lemma 1 m fcpa +i = m^pa mod p a+1 for every prime number p and k,a £ Z>o- 
Proof: The Laurent polynomial W(xi, . . . , x n ) has coefficients in Z. Therefore 

W{xi, . . . , x n ) k P a+1 = W{x{, xlfv" mod p a+1 Z[xf\ x* 1 ] . 
The lemma follows by taking constant terms. I 

Theorems 1.1, 1.2, 1.3 in l] together with the above lemma immediately 
yield the following integrality result for series and product expansions: 

Corollary 1 For z£R, z> C 2 

z-'exp^rj = z-' + ^AkZ-"- 1 = z-'l[(l-z- k )-^ (25) 
with Ak, bk £ Z for all k > 1. 



k 

,k>l I k>l k>l 



Remark. In 011] the result of Corollary ^ is used to interpret zQ(z) as the 
Artin-Mazur zeta function of a dynamical system, provided the integers bk are 
not negative. We have not yet found such a dynamical system within the present 
framework. 

For N £ N the moments of the measure dV/v(r) are, by definition, 



m W := / r k dV N (r)=N- n J2 w ( x ) k 
These are the same numbers as in l[17[). 
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Proposition 4 With the above notations we have 

m ( k N) > m k > for all N, k , 

nrifc = m fc if N > fcmax a! bg9maxi< 3 '< n \a,j - bj\ . 

Proof: N n (m^ — rrife) is the sum of the coefhcients of all non-constant mono- 
mials in the Laurent polynomial W{x\, . . . , x n ) k with exponents divisible by N. 
Since all coefficients of W{x\, . . . , x n ) are positive, this shows mjr^ > > 0. 

Assume N > fcmax a! bea. i<j<n \&j — bj\. Then all exponents in the mono- 
mials of the Laurent polynomial W(x\, . . . , x n ) k are > — N and < N. So only 
the exponent of the constant term is divisible by N. Therefore mjj. = rrife. ■ 

Note the natural interpretation (and proof) of this proposition in terms of 
closed paths on the graph L^r: closed paths on Fjv which are too short are in 
fact projections of closed paths on T. 



Corollary 2 For N > £ max max la, — b«\ and \z\ > C : 

a,be»l<?<n J J 



Q(z) ■ \B N {z)\ 



N~ 



exp 



This not only gives an estimate for the rate of convergence of 0) with respect 
to the usual absolute value on C, but it also yields the following congruence of 
power series in z~ x : 

B N (z)- N ~ n = 



~ 1 + J2 A ^~ k ~ 1 modz-'" 1 



fc>i 



with A 



k as in 



Remark about the relation with the large complex structure limit. 

Since the measure dV(r) is the push forward of the measure dt\ dti . . . dt n on 
W 1 by the function |I?(t)| 2 , one can rewrite (|13() as 

^ _ 1 f 1 dx\ dx 2 dx n 

(2ni) n J „ z- W(xi,...,x n ) x 1 x 2 " ' x n 

for z G C \ X. From this (and the residue theorem) one sees that H(z) is a 
period of some differential form of degree n — 1 along some (n — 1)- cycle on the 
hypersurface in (C*)" given by the equation W(x\, . . . , x n ) = z. As z varies we 
get a 1-parameter family of hypersurfaces. The function H(z) is a solution of 
the Picard-Fuchs differential equation associated with (that (n ~ l)-form on) 
this family of hypersurfaces. The Picard-Fuchs equation is equivalent with a 
recurrence relation for the coefficients in the power series expansion l|24|l of 
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H(z) near z — oo. All this is standard knowledge about Calabi-Yau varieties 
near the large complex structure limit and there is an equally standard algorithm 
to derive from the Picard-Fuchs differential equation enumerative information 
about numbers of instantons (or rational curves); see for instance |2(JI 1131 ITS) . 

On the other hand, we have the enumerative data nnjr^ of the present paper. 
In the limit for N — > oo these yield the moments and, hence, the Picard- 
Fuchs differential equation and eventually the instanton numbers. 



6 Examples 

6.1 n = l 

Mahler measures of one variable polynomials have a long history with many 
interesting results; see the introductory sections of 0] El- We limit our 
discussion to one example, without a claim of new results. This simple, yet 
non-trivial, example has n = 1, 21 = { — 1,1} C Z, c_i = Ci = 1 and hence 

|X>(t)| 2 = 2 + 2 cos(47rf) , W(x) = (x + x- 1 ) 2 . 

The moments are 

rrife = constant term of (x + x~ 1 ) 2k = 



2k 
k 



and hence by $Zfy: for z e R, z > 4, 



H(z) = £ ' 21 



k>0 



z -k-l 



k J V^4)' 



Q W = ex P ^y 7 ^=j=i(,-2-V^)) 
Applying Formula 117(1 to the present example we find 

"T= E ( 2 - 

j=k mod N 

which nicely illustrates Proposition 

Setting z = 2 + u + aT 1 one finds for the polynomials Bpf(z) defined in (J3J: 

BtfC*) = II (z-2-x 2 -x- 2 ) = u~ N Y[ (u-x 2 )(u-x- 2 ) 
= u N + u - N -2 

= -2 + 2^E( 2j 4^ -2)^-W. 
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So, Bn(z) is up to some shift and normalization the N-th Cebysev polynomial. 

The above computation also shows Bn(z) = Q(z) N + Q(z)~ N — 2 and thus, 
in agreement with (JIJ, 

lim B A r(z)- Ar_1 = Q(z) . 

N—*-oo 

For actual computation of Bn(z) in case z G Z one can use the generating 
series identity: 

rpN f T \ 

EM^ = -log l-(.-4)--_ . 

N>1 v \ I J 

For z = 6 one finds (using PARI) 

B N [Q)T N = 2T + 12T 2 + 50T 3 + 192T 4 + 722T 5 + 2700T 6 + 10082T 7 

+37632T 8 + 140450T 9 + 524172T 10 + 1956242T 11 
+7300800T 12 + 27246962T 13 + 101687052T 14 
+379501250T 15 + 1416317952T 16 + 5285770562T 17 + . . . 

For primes p in the displayed range the number i? p _i(6) is divisible by p 2 for 
p = ±1 mod 12 and is not divisible by p for p = ±5 mod 12 and is exactly 
divisible by p if p = 2,3. We also checked 5 2 | B 24 {6) and 7 2 | S 48 (6). This 
agrees with the number of solutions of the equation u + u^ 1 = 4 in ¥ p and F p 2 . 

If z € Z, z > 4, then Q(z) = 5(2; — 2 — \J z{z — 4)) is a unit in the real 
quadratic field Q(yJ z(z — 4)). According to Dirichlet's class number formula it 
relates to the L-function of this real quadratic field: 

log(QW) = ^L(l, X ) 

where D, h, \ are the discriminant, class number, character, respectively, of the 
real quadratic field Q(y/ z(z — 4)) (see e.g. The relations between Mahler 
measures and values of L-functions, which have been observed for some curves, 
are perfect analogues of the above class number formula (see [T5]V 

6.2 The honeycomb pattern. 

For a nice two-dimensional example we take 21= {(1,0), (0,1), (—1,-1)} C Z 2 , 
c (i.o) = c (o.i) = c (-i,-i) = 1 an d hence 

\V(t 1 ,t 2 )\ 2 = 3 + 2cos(27r(ii - t 2 j) + 2 cos(27r(2ii + t 2 )) + 2 cos(27r(ii + 2t 2 )), 
W(xi,x 2 ) = (xi + x 2 + x^ 1 x^ 2 1 ){x^ 1 + x^ 1 + X\X 2 ) 

= x\x 2 x + x\x 2 + x^ x x 2 + x\x\ + Xi 2 x 2 x + Xi 1 x 2 2 + 3 . 
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As a basis for the lattice A we take (2, 1) and (—1, —2). This leads to coordinates 
u\ = x\x 2 and u 2 = x% on the torus U A . In these coordinates the function 
W reads 

W(ui,u 2 ) = ui + Ui 1 + u 2 + «2 + u^ 1 u 2 + uxu^ 1 + 3 (26) 
= (iti + u 2 + 1)(% X + U2 1 + 1) . 

Figure |21 shows a piece of the graph T. Figure shows some level sets of the 
function \T>(t 1 ,t 2 )\ 2 ■ The dual lattice A v is spanned by (1,0) and (§,§). The 
maximum of the function \D{ti, t 2 )\ 2 equals 9 and is attained at the points of A v . 
The minimum of the function |X? (i x , ^2 ) | 2 equals and is attained at the points 
of (0, — |) + A v and 0) + A v . There are saddle points with critical value 1 at 
(— g, ^)+A v , (i, — i)+A v and (g , i)+A v . In terms of the coordinates u\, u 2 the 
maximum lies at (u\, u 2 ) = (1, 1), the minima at (e 2 ™ /3 , e 4 ™ /3 ), (e 4 ™ /3 , e 2 ™ /3 ) 
and the saddle points at (1,-1), (—1,1), (—1,-1). The algebraic geometry of 
this example concerns the 1-parameter family of elliptic curves with equation 
z — W(ui, u 2 ) = 0. In homogeneous coordinates (Uq : U\ : U 2 ) on the projective 
plane P 2 , with u\ = UiU^ 1 , u 2 = U 2 Uq 1 , this becomes a homogeneous equation 
of degree 3: 

(UoUi + U U 2 + C/i[/ 2 )(C/ + U! + U 2 ) - zU a UiU2 = . (27) 

Beauville 2 showed that there are exactly six semi-stable families of elliptic 
curves over P 1 with four singular fibres. The pencil l|27|l is one of these six. It 
has singular fibres at z = 0, 1, 9, 00 with Kodaira types I 2 , J3, I\, I§, respectively. 
Note that the first three match the critical points and levels in the diffraction 
pattern. After blowing up the points (1, 0, 0), (0, 1, 0), (0, 0, 1) of P 2 one gets 
the DelPezzo surface dP3. The elliptic pencil (|27|l naturally lives on dP3. It 
has six base points, corresponding to six sections of the pencil. Since the base 
points have a zero coordinate, these sections do not intersect the real torus U A . 
Equations (|26(l and l|27|l also appear in the literature in connection with the 
string theory of dP3. 

Formula Q23JI and some manipulations of binomials give the moments: 




These numbers satisfy the recurrence relation (see ^1 Table 7) 

(k + l) 2 m fe+1 = (10/c 2 + 10k + 3)m k - 9/c 2 m fe _i . (28) 

We refer to J5| Example C and to [TS] Example j}6 for relations of these numbers 
to modular forms and instanton counts. Golyshev |1(J| can derive the recurrence 
1281) from the quantum cohomology of dP3. The numerical evidence for the 
relation (|21ll between Mahler measure and L-function in this example is given 
in gj Table 2. 
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With Formulas 117|) and (|26[) one easily calculates 



m 



(JV) 



E 

ii=i2 mod N,ji=j2 mod N 



k 

ii 

AN) 



k — i\ 
h 



k 



k- i 2 
h 



Note that these formulas confirm = rrifc for k < N. 

Equation l|27[l is clearly invariant under permutations of Uq,U\ 1 U 2 . There- 
fore the diffraction pattern has this ^-symmetry too. Since only the critical 
points have a non-trivial stabilizer in S3 the multiplicities mult tv (r) in this ex- 
ample satisfy 



mult 7v (9) 
mult 7v (0) 
mult 7v ( 1 ) 
mult tv (r) 



1 

2 

3 mod 6 
mod 6 



VN 

if 3|JV 
if 2|JV 
if r + 0, 1, 9, 



(29) 



ViV. 



We have computed the numbers multTv(r) for some values of N. We found for 
instance 

B 6 (z) = z 2 (z - l) 15 (z - 3) 6 (z - 4) 6 (z - 7f {z - 9) . 

We computed W(ui,u 2 ) = («i + u 2 + ^(uV 1 + u^ 1 + 1) for ui,u 2 € F7: the 
(i,j)-entry of the following 6 x 6-matrix is W(i,j) mod 7: 



3 

3 4 

4 

1 

1 4 



3 1 

1 1 

1 4 1 

3 4 1 

4 1 



11111 

This yields the following count of points over F7: 



mod 7 







1 
15 



3 4 5 6 
6 6 



BUG (¥* 7 f I W{0 =zinF 7 } : 
Thus we see that the inequality in (|20() can be strict: 

u 7 (B 6 (53)) = 12 > 6 = m e (F 7 ) 2 I = 53 in F 7 } . 



(30) 
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